Renormalization of resonant tunneling in MOSFETs 
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We study tunneling between a localized defect state and a conduction band in the presence of strong 
electron-electron and electron-phonons interactions. We derive the tunneling rate as a function of the 
position of the defect energy level relative to the Fermi energy of conduction electrons. We argue that 
our results can explain the large tunneling timescales observed in experiments on random telegraph 
signals in Si metal-oxide-semiconductor field effect transistors. 



PACS Numbers: 73.40. Qv, 63.20.Mt, 85.30.Tv 

It has been long realized that two-state systems can 
have profound effects on micro- and nano-scale electronic 
devices, leading to random switching or apparent low- 
frequency noise in transport. In the well studied case 
of Random Telegraph Signals (RTS) in Si Metal-Oxide- 
Semiconductor Field Effect Transistors (MOSFET), elec- 
trostatic measurements reveal that the switching is most 
likely caused by the fluctuating charge of interface de- 
fects [1]. The charge fluctuations on the defects lead to 
fluctuations in the scattering potential experienced by 
the carriers. In small devices this leads to bi- or multi- 
valued fluctuating conductance [1], while in devices with 
a large number of defects it results in l//-like frequency 
dependent noise [2]. 

At low temperatures, the charge fluctuation of the de- 
fects are mediated by quantum mechanical electron tun- 
neling between the localized defect states and the itin- 
erant states. The signature of the quantum-mechanical 
tunneling is the temperature-independent switching rate 
at low temperatures [3]. In this problem, we must also 
consider the Coulomb interaction, which in fact makes 
RTS observable. Tunneling in the presence of strong 
electron correlations has been of active interest as it ex- 
hibits a number of peculiar features. One of these fea- 
tures is the Fermi-edge singularity in the resonant im- 
purity tunneling. The singular dependence of the tunnel 
rate on the energy of the impurity electron relative to the 
Fermi energy of conduction electrons was first suggested 
by Matveev and Larkin [4], whose theory was based on 
the X-ray emission/adsorbtion problem. The latter was 
extensively studied in numerous works, see Ref. [5] and 
references therein, which have shown that the transition 
rate is proportional to , where E is excess energy rel- 
ative to the absorption/emission threshold. Nozieres and 
De Dominicis [6] have demonstrated that the power a 
is related to the scattering phase shift of the conduction 
electrons at the core hole potential and it has two es- 
sential contributions: a negative one due to exciton-like 
physics [5], i.e. an attraction between the core hole and 
excited electron, and a positive one due to the orthog- 
onality catastrophe [7], i.e. an adjustment of the Fermi 
sea to the appearance of the core hole potential. For a 
single channel model, i.e., in case of spinless Fermi sea 



electrons interacting with point-like core hole potential, 
Nozieres and De Dominicis obtained that 
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where XiWC^ivV) is the scattering phase shift, v is the 
density of states of the conduction electrons and V is 
the Coulomb interaction strength. In the case of a local 
defect tunnel- coupled to a conduction channel, the same 
model applies after making an association between the 
X-ray electrical field coupling strength Ed and the tun- 
neling matrix element A, and defining the energy E as 
the difference between the localized level position and the 
Fermi level. Then, one finds that the tunnel rate between 
the impurity and the Fermi sea is [4]: 
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where 70 = 27ri'A^, O is a step function, Tg is the cutoff 
of the order of the conduction electron Fermi energy and 
a is given by Eq. (1). 

The singular behavior predicted by Eq. (2) has been 
observed experimentally in the resonant impurity tun- 
neling in Si MOSFETs [3], with a negative exponent a. 
There is, however, a serious discrepancy between the sim- 
ple theory prediction, Eq. (2), and experiments [3,8] that 
has to do with the prefactor 79. In particular, in Si MOS- 
FETs the tunneling is found to occur on a rather slow 
timescale, typically from milliseconds to seconds. 

(a) Met. Ins. Semi. (b) Ins. Semi. 
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FIG. 1. (a) Schematics of MOSFETs; (b) Energy diagram 
for the impurity level and the conducting electrons. Coupling 
to optical phonons shifts the bare position of the level. 

Given such tunneling timescales one might have 
guessed that the slow tunneling rate is due to a large 
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distance, on the order of 15-20 A, between the impuri- 
ties and the insulator-semiconductor interface. However, 
the distance between the resonant impurity and the 2- 
Dimensional Electron Gas (2DEG) in the MOSFET con- 
ducting channel can be determined via a quite simple ex- 
perimental procedure based on the monitoring of the de- 
fect average occupancy as a hmction of the gate voltage, 
which reveals that the separations between the impuri- 
ties and the 2DEG is typically of the order of 1-2 A [11]. 
A naive estimate suggests that for such width of the tun- 
nel barrier (of height of several eV) the impurity electron 
dwell time must be in the range of pico -to nanoseconds, 
which is by many orders of magnitude smaller than the 
observed dwell time. In the present work we argue that 
this apparent contradiction can be resolved if one takes 
into account strong electron-optical phonon coupling in 
the Si MOSFET insulating layer, which drastically slows 
down the defect charge dynamics. To this end, we solve 
the problem that simultaneously includes the effects of 
electron-electron and electron-phonon interactions. We 
find that at energies E small compared to the optical 
phonon frequency one recovers the singular behavior of 
Ref . [4] but with a renormalized "bare" tunneling rate 70 , 
while at sufficiently large energies (for E > Ep, where Ep 
is the polaronic shift defined below) the renormalization 
vanishes. We make numerical estimates which indicate 
that the renormalization effect can be very strong in a Si 
MOSFET. 

A schematic diagram of MOSFET is presented in Fig. 
1(a). The impurity is located near the interface in the 
insulating layer. The Hamiltonian of the system includes 
the electronic part Hq and the electron-phonon interac- 
tion Hph- The electronic part is 

Ho = Y^ Ekclck + EodU + VdUY^ clck' 

k k,k' 

+AJ2{cid + d^Ck). (3) 

k 

Here c\,{ck) and rft(rf) arc crcation(annihilation) opera- 
tors of spinless electrons in the Fermi sea and at the 
impurity respectively (we will discuss the role of spin 
degeneracy later on), the third term describes Coulomb 
coupling of the impurity electron and the Fermi sea elec- 
trons, and the last term represents tunneling between the 
impurity and conduction electrons. We assume that the 
tunneling matrix element A is independent of the single 
particle states in the Fermi sea and that the impurity 
potential is point-like. We also assume that the density 
of states (v) of Fermi sea electrons is constant and set 
the chemical potential (/x) in the Fermi sea to zero in the 
following calculations. 

In general, electron-phonon interaction in bulk semi- 
conductors is insignificant. However, in the insulating 
oxide layers of MOSFETs, the polar structure of the ma- 
terial favors strong electron-phonon coupling. We model 



this local electron-phonon coupling by a Holstein-type 
Hamiltonian 

-ffph = g d'' d {a' + a) + ujQa^a , (4) 

where o^(a) are creation(annihilation) operators for a lo- 
cal optical phonon of frequency ujq coupled to the occu- 
pation number at the impurity. The coupling constant g 
will be estimated later [6]. 

We now evaluate the defect electon tunnelnig rate at 
T = by calculating the matrix clement (0|(!^'*"^|0), 
which is the probability amplitude for the impurity state 
to remain empty. Here H = Hq + Hph and the "ground" 
state |0) corresponds to the trap being empty, and the 
phonon and the conduction electrons in the ground state, 
i.e., |0) = |0)d (E> |0)ph fSi |0)ci. We will carry out the 
calculation in Euclidean time by setting it = (3. In what 
follows h and electron charge e are set to unity unless 
stated otherwise. 

We expand Z — (0|e^''^lO) in terms of the tunneling 
term in the Hamiltonian H as 

^ = E / d^''[(i]me-^"'HT{M...HT{P^)\Q), (5) 

N •' 

where we have introduced a shorthand notation 
jTd^^m = !^dP2N!t" dP2N-i-!tdPi. r is the 
time ordering operator, and N = 0,1,... . In 
Eq. (5) HtW) = exp {(}H')Ht exp {-I3H'), where Ht = 
A Y^};{c\d-\-d'^ Ck) and H' contains all the remaining terms 
of the Hamiltonian H. Note that H' is diagonal in d'^d 
and can be decomposed into two commuting parts, i?ei 
and ffph, that involve conduction electrons and phonons 
respectively. 

The conduction electrons and optical phonons can be 
integrated out by realizing that a term with given A'' in 
the sum in Eq. (5) corresponds to trap site being occupied 
for periods of time /32i-i < (3 < Pn, i < N , and being 
empty otherwise, see Fig. 2(a). One can formally write 
that term in the sum in Eq. (5) as 

N 

A^^(0|Tn E 4(/32,)cfc.,_.(/32,-i)x 

j=l k-2j,k-2j-i 

^-/>e.[n.(/3')]<i/3'|o^^^ ^0|^^-/>,U".(/3')]c^/3'|0)p^. (6) 

Here the explicitly time dependent Hamiltonians 

Hci[nd{Py\ and H-pi\nd{P)] correspond to the first three 
terms in Eq. (3) and to Eq. (4) respectively, with the re- 
placement d'^d — > nd(/3), where nd{(3) is the occupation 
number of the trap as a function of time, 

N 

naiP) = E - " ©(/^ " /?2i)] ■ (7) 

The phonon-dependent matrix in Eq. (6) can be evalu- 
ated by a standard method of either cummulant expan- 
sion or by path integral techniques as a ground -to ground 
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state amplitude of a harmonic oscillator subjected to an 
external force nd,{(3). By using Eq. (7) this matrix ele- 
ment is 



, (8) 



j>k 



where we have defined the polaronic shift Ep = .g^/wo 
and G = Ep/ujQ. The first exponent in Eq. (8) is due to 
the diagonal terms (j = k) in the sum over j, k. 

The electronic matrix element in Eq. (6), which is es- 
sentially the iV-particle propagator for the conduction 
electrons, can be evaluated in the large /3 limit by the 
technique developed by Nozieres and De Dominies [6], 
and by Anderson and Yuval [12]. The technique relies 
on the flat density of states assumption for the Fermi sea 
electrons (which is the case for the 2DEG in the con- 
duction channel) and is based on the theory of singular 
integral equations. It allows one to obtain an asymptoti- 
cally exact expression for the A^-particle Green's function 
in Eq. (6) in the limit \P2j — P2j~i\ ^ ^e- By adjusting 
the result of Refs. [6,12] for the electronic Green's func- 
tion in Eq. (6), and substituting Eq. (8) in Eqs. (5,6) 
we obtain 



N 



7oe ^ 
27rrp ' 



exp^(-l)-''-^-[Xln(^^^) - Ge-^o{P,"^^^] . (9) 

In Eq. (9) we have defined Er = Eq — Ep, — , 
and K = a. Thus, due to the coupling to the phonon 
mode the resonant level shifts downwards as shown in 
Fig. 1(b). 



(a) 



P2J-1 P2j 



P2j-1 02j+2 







domain walls 










'i i U i i i i i i 







FIG. 2. (a) Typical "trajectory" for the trap occupation 
number as a function of time; (b) Spin domains in the equiv- 
alent Ising model 



Fig. 2(b). The spacing in imaginary time between suc- 
cessive spins equals Tg. This Ising model has a nearest 
neighbor coupling term — (l/2)[ln (27r/7ore) -I- G\siSi^\ 
{si = ±1), and a long range interaction —{1/2)[K [i — 
j)^^ + {GuqtD exp{—u)oTe\i — j\)\siSj. The renormalized 
position of the resonant level En plays the role of an 
external magnetic field. In the weak long-range inter- 
action regime relevant here, K < 1, the Ising system is 
in a paramagnetic state composed of ferromagnetic do- 
mains whose size is controlled by the short range interac- 
tion. A simple estimate for the mean separation between 
the domain walls at zero external field {Er = 0) yields 
I « exp (G')/7o ^ Te- The external field {Eji ^ 0) forces 
spins to align along its direction, thus shrinking the do- 
mains opposing the field relative to the ones aligned with 
the field. The ratio of the corresponding domain sizes 
can be estimated by evaluating the average magnetiza- 
tion (s) per spin. For E^ ^ 70 exp (— G/2) one gets 
1 — (s) ~ 7o exp {—G)/{TeE]^) and therefore the average 
size of the domains with spins along the field exceeds that 
of the domains of spins opposing the field at least by a 
factor of 77 = rei?|./[7o exp (— G)] ^ 1. (The long range 
terms will favor even more polarized state as they pro- 
vide additional ferromagnetic coupling.) As a result the 
interaction between domain walls from different domains 
is negligible provided loq, En ^ 70 exp [—Ep/uo]. 

The above arguments allow for a significant simplifica- 
tion of the transition matrix element Z in Eq. (9) in the 
limit of small joTe exp (— G). Retaining only those terms 
in the last exponent of Eq. (9) that couple f32j and f}2j-i 
we obtain that Z = 



E 

N 



70 
27rr„ 



N 



N 

/ d2^[/3] n</'(/32,-/32,-i) , (10) 



where (/)(/?) = exp [-ErI3 -Kin {fi/Te) + Ge-"'"^]^ q^^^. 
approximation is similar to the non-interacting "blip" ap- 
proximation used in the solution of the spin-boson prob- 
lem [14] and therefore the same method of Laplace trans- 
form can be used to evaluate the sum in Eq. (10). Upon 
Laplace transform Z = d(5e~^^ Z{[3), which converts 
theA^'s term in the sum into (Ae-^)2^(l/A)^+i[^(A)]^, 
where 0(A) — d(3e~^^(f){(3). Re-summation in 
Eq. (10) then yields Z = [A - (7o/27rTe)e-G0(A)]-i. The 
tunnel rate 7 can be associated with twice the imaginary 
part of the pole of Z. In the limit 7oTee~'^ ^ Er, the 
rate 7 of tunneling into the impurity is given by 



-G 



2-KTe 



Im / exp 



Te 



dp, (11) 



To proceed, we write Eq. (9) as the T — 1 parti- 
tion function of a one-dimensional ferromagnetic Ising 
model, where a spin up or down corresponds to an 
empty or occupied impurity respectively, as shown in 



In order to evaluate the imaginary part of the above in- 
tegral an analytical continuation must be used. We fol- 
low prescription of Ref. [15]. The procedure is straight- 
forward: at the saddle point /3o (given by equation 
—Er — {K/(3o) = Gexp (— wq/^o)) we deform the contour 
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of integration into the complex plane in the direction 
of steepest descent (along the +i direction in our case). 
As a result the integral picks up a convergent imaginary 
part. The numerical result is presented in Fig. 3. In 
the limiting cases of small and large En simple asymp- 
totic expressions for the tunnel rate into the impurity are 
possible: 



7in = T'o^2^Q( ^^^ lERTe^ , \Er\ « ivo I (12a) 



r(l +a) 



r(l + a) 



[iER-Ep)Ter, \Er\>Ep. (12b) 



In Eqs. (12) f(z) = (27r/z)i/2z^e-^ For Er < Ep the 
rate is suppressed by exp {— Ep / ujo) ■ This suppression 
occurs in case of elastic tunneling, when no real phonon 
modes are excited, and essentially represents the overlap 
of the phonon ground state wave function for occupied 
and empty defect states. [14]. The exponent a is due 
the X-ray singularity as in Eq. (2). For Er ^ Ep the 
exponential rcnormalization disappears due to opening 
of additional inelastic tunneling channels with excited 
phonons in the the final state. 

So far we have not included the effects on spin degen- 
eracy. However, based on the close analogy between our 
results Eqs. (11,12) and the X-ray absorption [6] prob- 
lem, we argue that the main effect of the spin degeneracy 
is to change the second term in the expression for the ex- 
ponent a in Eq. (1) by a factor of 2. 




FIG. 3. Dependence of the tunnel rate on the bias 



For the interface impurity states in a Si MOSFET 
the coupling constant turns out to be sufficiently strong. 
In order to estimate the order of magnitude of the po- 
laronic shift that enters Eqs. (12) {Ep = g'^/ujo), we 
take Frohlich coupling constants of an electron to an 
optical phonon mode with a wavevector k [13], gk = 



a(*d|e''=''|^'d)/(V"i/2|/j|)^ Here V is crystal volume and 
= 27re^wo(e^^ — e^^), where eoo and eo are high 
frequency and static dielectric constants, e is electron 
charge and ^'^ is the wavcfunction of the impurity 
electron. Assuming that the latter is a hydrogen-like 
ground state with effective Bohr radius Cd, one finds that 
(*d|e*'"'|*d) = [1 + (adfc)^/4]"^. Using the equivalence 
between electron-phonon coupling in our single-mode 
model and Frohlich's multi-mode Hamiltonian [10] , that 
is g"^ = J2k9k^ "^"2 obtain that Ep = 5e^/(16ad)(e^-^ — 
e,^^). For a deep impurity level in Si02, ~ lA, Coo — 2, 
eo — 4, we get Ep ~ leV. The frequency of the opti- 
cal phonons near the Si02 — Si interface is of the order 
60 meV, which yields the renormalization of the tunnel 
rate in Eq. (12a) by a factor of exp (—16) ~ 10"''. While 
the detailed quantitative theory warrants further study, 
we believe that our model gives a qualitative explanation 
for the tunneling slowdown at Si — SiOg interfaces. 
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